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(57) Abstract 



A method for determining a physical or 
atavistic of a liquid fat which a drop of the BqnJd is 
on the surface of a crystal of a is— ti , crj 
The liquid is evap o rate d from said sorface vhfle a 
put of the crystal migobalancn hi ineasnsedL Ileal 
is amuyaed to determine me dtamiteiisoe of the Bcjafci 
inay be. foreramia^ 
ity. Apparatus for carrying out the method is 
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HETH00 AND APPARATUS FOR DETERMINING A PHYSICAL OR CHEMICAL CHARACTERISTIC OF A 
LIQUID 

The present invention relates to the determination of the characteristics of a liquid in 
particular by using a quartz crystal microbalance (QCM). The characteristics may be, 
for instance, the chemical composition or the viscosity of the liquid 

5 

The QCM was first used as an accurate weight measurement device for measuring 
die weight of thin film metals. The microbalance works by applying an oscillating 
electric field across a quartz crystaL The field causes a shear oscillation in the crystal 
known as the converse piezoelectric effect and the crystal oscillates at a stable 

10 resonant frequency. When matter is deposited cm die crystal the effective wavelength 
of the oscillation is increased and the frequency is thus reduced Therefore, the 
resonant frequency of the crystal is sensitive to matter deposited on its surface or 
which is coupled to the surface by viscoelasric means. In addition when connected to 
a network analyser, admittance, quality factor Q, phase and radio frequency (rf) 

IS voltage may be measured The effect of various interfecial physical phenomena on 
these quantities has led to die use of the QCM in many analytical applications 
including adsorption of gases, the solid-liquid function, electrochemistry, thin films, 
polymers and the study of biological reagents at the crystal surface. 

20 The me of the QCM for lneamcm ent of fluids has been explored In c omp a ris on 

with a so lid mass, which is entirely mechanically coupled to the quartz surface due to 
its rigidity, die reduction in resonant frequency in a fluid is due to an effectively 
rigidly-coupled surface layer. Beyond this layer, the crystal is unaffected by the fluid 
other than by the small effect of the pressure exerted by the weight of the fluid 

25 

As a consequence, analysis of fluids has concentrated cm measurement of mass 
change in a stable liquid layer on the surface of die miczobalance. Such a fluid is in 
equilibrium and can provide a uniform coating across the entire crystal surface. 
Changes in the comp o siti on of the fluid can therefore be determined by measuring 
30 the change in mass per unit surface area. 
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According to an aspect of this mventico there is provided method for determining a 
physical or cheniical characteristic of a liquid, the method comprising: 
depositing said liquid on a surface of a crystal of a quartz microbalance to form a 
droplet on said surface; 
5 evaporating said liquid from said surface; 

measuring a signal output of the crystal microbalance; and 

analysing of the signal output to determine said characteristic of said liquid. 

The droplet on the surface of the crystal may be described as a sessile body of liquid, 
10 indicating that it simply sits on the crystal surface. 

The characteristic of the liquid may be, for instance, some aspect of the chemical 
composition of the fluid or its viscosity. 

IS Evaporation of die liquid drop produces changes in a number of measurable physical 
characteri stics of the crystal any of which may be m easured as a signal output from 
the crystal. The measured values are characteristic of the chemical composition of the 
liquid. 

20- Thadnrnktimfinbtylinavolia^ 

Preferably, the output signal is dependant on the evaporation rate of liquid from the 
crystal surface. 

25 Preferably the change in resonant frequency (4£)is measured. Optionally, 4^ can be 
measured by incbding the crystal microbalance as a cornponcnt in an oscillator 
circuit. 

Optionally, any or all of the following physical chai aul e i isn cs, namely, admittance, 
30 quality factor (QJ, phase or radio frequency (rf) voltage may be measured. 
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Preferably the signal ismeasured as afunction of time. Optionally the 
are taken at least three times per second. 

Preferably the measurable physical characteristics may be measured using a network 



Preferably the liquid may be composed of a single uiitaio^ chemical or a mixture of 
unknown chemicals. Optionally, the liquid may be an organic liquid. A further 
option is that the liquid mixture cootaii^ 



10 

Preferably, the signal output is analysed by comparmg values derived f^ 

output wife known values contained in a database. Such a comparison allows the 

unknown liquid to be identified. 

15 The present invention also provides an apparatus for d et e rminin g a physical or 
chemical cliaracteristic of a liquid, the apparatus coinpri^ 
a quartz oystal mtcrohalaoce; 
means for depositing said liquid on the surface of said 
droplet on said surface; 

for measuring the signal output of the crystal mkrobalance; and 



For many volatile liquids the observed frequency response is qualitatively 
char acteristic of the specific liquid used such that a noo-deterministic recognition 
method can be employed in order to identify the Kqu^ The origin behind such 
characteristics is believed to be the variety interaction between convection and 
conduction processes, and the surface tension of the liquid. Such liquids can be 
termed unstable-interface liquids. However, tor other, staMe-haetfece liquids, the 
responses are generic, differing only in the severity of the response curve. The 
response curve is paraineterised by the crystal se^ 
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and ft is the angk subtended from the crystal centre across its surface This sensitivity 
is formally expressed as a series of Bessd functions but is more often approximated 
as Gaussian, as in equation (1): 

5 5(r,#)-S t ex^.W)^j (1) 

where /T controls the sensitivity dependence on radius r and the radius of the QCM 
electrode is r+ it is possible to relate /jf(t) to elapsed time At as provided in the linear 
expression of equation (2): 

where Afma is the maximum change of frequency observed, r(U) is the radius of the 
drop at the moment it is deposited and v r is the retreat speed. By dividing through by 
IS -4f7(]-exp-fi) and taking logarithms: 




25 where v, is defined in terms of area change per unit time, m s" . 
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Preferably a computer is used to store a database of known values. Preferably a 
computer is used to compare the known and unknown values and select the best fit 
from the known signals to determine die identity of the unknown chemical. 

5 Preferably the crystal is driven at its r eso nan t frequency or harmonics thereof by an 
Alternating Current supply. Preferably the crystal nncrobabnee is constructed from 
quartz crystal. Preferably the quartz crystal is unpolished. Optionally the quartz 
crystal has a diameter of between 5mm and 15mm. 

10 Preferably the quarts crystal microbalance is attached to a network analyser by 
means of an electrode on both its upper and lower surfaces. 

Preferably die rate of evaporation is controllable. Preferably the rate of evaporation is 
controlled by controlling die temperature of the crystal surface. Optionally, the rate 
15 of evaporation is controllable by controlling the pressure at the crystal surface. 

Preferably, die volume of liquid deposited on the surface can be controlled 

The method of the present invention can be used to detennine the viscosity of die 
20 liquid droplet. From die Sauerbrey equation we know that the change in oscillating 
frequency of a Quartz crystal microbalance transducer is related to the mass loading 
that crystals surface (equation (1> 

y -2/. 2 Am 

(1) 

30 

A/= frequency change, / Q = crystal r eso nan t frequency, A*r= mass change, A= 
electrode area, = shear modulus of quartz and /y= density of quartz. 
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When a dioplel is in contact with an oscillating sirface, Ibc owllMions passing from 
the surface into the fluid droplet decay according to an exponential law. The 
penetration depth 8 is a distance normal to the crystal surfecc over which the 
amplitude of the oscillations reduces by e times. This penetration depth is dependent 
S on the properties ofthe fluid forming the droplet as seen m equation (2). 



40 



15 where:- ijm" absolute fluid viscosity, and /»= fluid density. 



The Sauerbrey equation assumes a solid mass loading of the crystal where the whole 
of the mass oscillates with the crystal frequency. As explained above, when the 
loading is fluid, the whole volume of the droplet will not oscillate with die crystal as 
20 the amplitude decays through the fluid. The volume of fluid oscillating at the crystal 
frequency is equal to a volume of A*tit2. This means it can be assumed that a layer 
forms on the electrode, which acts as a solid mass and is known as the rigidly 
c«mkd layer. This layer is the Am res 

the crystal iJr"itinm. The mass of this layer b related to the density and volume of 
25 the dropfct as shown in equation (3). 



30 Am=p,A^ (3) 

Substituting equations (2) and (3) into equation (1) gives an equation relating 
35 frequency change to viscosity. 



This rearranges for viscosity as shown in equation (5) 
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The invention will now be described by way of example with reference to the 
accompanying drawings of which: 

Figure 1 shows die equipment and experimental set up used for measuring 
10 liquid characteristics on evaporation; 

Figure 2 shows a set of graphs for plotting the change in resonant frequency 
of the crystal against time for a range of liquids on evaporation; 

Figure 3 shows a set of graphs plotting the change in resonant frequency 
against time for a liquid using different types of crystal; 
15 Figure 4 show s a schematic diagram of a device for determining die content 

of a liquid; 

Figure 5 is a graph showing the change in resonant frequency against time for 
butan-l-ol; 

Figure 6 shows response curves for a range of alcohols; and 
20 re prtrt 7 Amw tfce results of viscosity mcaauiciucnts on TBP/OK mixtures. 



Referring to Figure 1, a quartz crystal was used and operated at a resonart treqiiency 
of 10MHz by frequency generator 5. The quartz was unpolished, with total diameter 
ff f 8mm «nrf « nrwer electrode of ap proximately 4mm diameter on each free. The 

25 crystals were connected to a network analyser 7 using grounded coaxial leads to 
mammae stray capacitance effects and external interference. The crystal surfaces 
were orientated in the horizontal plane. The network analyser 7 used in this work 
was a Hewlett Packard 8753C and was interfaced to a personal computer 9 via Lab 
View, a data acquisition software package. The network analyser 7 was set up to 

30 record the change in frequency Af from the resonant frequency/ at a rate of 3 



In performing the experiment a droplet of alcohol 3 was applied to upper surface of 
the crystal 2 using a syringe (not shown). The droplet was observed to completely 
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cover the upper suite 2 of the crystal and a men droplet vohinm of 1mm was 
^corded. Data were recorded from a time immediately prior to the droplet being 
deposited on the upper surface of tte.crystal2toiinmedn^ 



selection of volatile alcohols were used m the expermients to study the effect of 
molecular chain length on over time. 



•Hie above experiment examined methanol, ethanol, propan-2-ol, butan-l-ol and 
10 pentan-l-oL 



Figure 2 shows the change in resonant fw,uency of Ai(Hz) agamst 
the above alcohols. 2 The plots (a) to (e) show the results for me&an^ 
propan-2-ol, butan-l-ol and pentanrl-ol, respectively. 

m general, the response of the oscfflatmg quartz crystu to 

cvaponrtingali^ Ahhougha 
srmnar general response b eaJriMted by aJI grve alcohol 



Fnnne 3 stows the m resonant f««^ 

oq three separate crystals of the same type. 



AJomugh uiere are changes m tto 
25 volume, which resuh m changes m tto le^ 

characteristic shape of resiwnrefarepwduced. This was evident for the other alcohol 



In Figure 4 a device 31*shownwhichcanbeusedto determine the cluneal 
30 contents of a sample. A senii»te 32 b placed on tto up 

39 by means of a syringe or other accurate nweuring device. Use of such an 
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accurate measuring device ensures thai a substantially constant volume of sample 32 
can be analysed. Temperature control dements 35 are used in conjunction with a 
thermometer (not shown) in order tp control the temperature inside the evaporation 
chamb er 33. T his in turn controls the rate of evaporation for the sample. As in 
5 example 1, the crystal is driven by an ac signal generator 49 at the resonant frequency 
of the crystal (or harmonics thereof). 

On evaporation, the change in frequency Afis sampled by a network analyser 41 
which is set to sample 4/~3 times per second. Sampling rate can be easily increased 
10 lo improve accuracy. The data obtained during sampling is then stored in a memory 
chip 43. Otoceallofthe data has been collected, the stored data is downloaded onto a 
computer where analysis of the data is undertaken. 

Figure S shows die change in resonant frequency Af (Hz) against time (t) for butan- 1 - 
IS ol with a least squares fit to the region where the evaporation rate is c on sta nt . This 
linear region of the graph is due to the increase in evaporation rate being 
co m pensa te d forbya decrease in mass on the crystal surface. After die linear section 
the change in mass dominates Figure 5 has a characteristic shape similar to that of a 
ienni luncaoo. 

Data analysis can be performed by using equation (3) set out above. The value of the 
constant fi has been taken to be 2.00±0.05 and mrirnendrnt of A which is consistent 
with other studies of radial sensitivity of quartz crystals. Plots of 
v=- h+SfQX} - exp- P) I V— - 1) against 4r for the series of alcohols studied are 
25 presented in Figure 6. The alcohols are a) methanol, b) ethanol, c) propan-2-ol, d) 
butan- l-ol and e) pentan-l-ol. 

. The data of Figure 6 are 1^ q^a^^ of 1^ caa^ iMm M mall 

cases, bodi the earfy data (t<l OOsX and that just prior to the crystal reaching #~0 Hz, 
30 exhibit variations fiom the central linearity winch is characteristic of radial 
sensitivity effects currently beyond the model used Indeed, for methanol the linear 
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feature between these 

this alcohol. 



non-linear aspects is very short due to the rapid evaporation of 



The region of each data set presented in Figure 6 has been fitted using a linear 
5 regression in order to extract the ndial retreat speed * These results are presented 
in Table 1: 



Alcohol 


Methanol 


Etfaanol . 


Propanol 


Butanol 




Region of fit 


95-105S 


96-99s 


6040s 


100-2005 


300-5008 


Degrees of 


29 


8 


59 


299 


599 


freedom ^ 












R 2 


0.94 


0.964 


0.994 


0.988 


0.953 


Slope flO* a 1 


1.16*0.05 


9.68±0.66 


?77H).02 


0.465*0.003 


0231*0.002 




-1.71±0.05 


-9.77*0.64 


-1.99±0.02 


-1.082*0.005 


-1353*0.009 


*/ir 7 mV 


436*840 


3841*239 


8.72±0.08 


1.83*0.01 


0.907*0.008 
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TABLE. 1. Results from the linear regression analysis of die data in Figure 6. 

In summary, a quantitative method of fluid identification is possible, in addition to 
the respoi« recognition method, vio the retina 
advanta geous in those cases where die frequency 
similar, for example those of butanol and pentanol in the date provided here. 



In cases where the liquid is unknown, the date availabte from 
change in resonant fiequency with evaporation rate can be correlated to a library of 
known values. The measurements can be subjected to signm ptocessmg, for mstaice, 
20 by Fast Fourier transformation. 

A further example of use of the above technique is tor the analysis of samples of die 
extractent tri-butyl phosphate (TBP). Despite being chenricaUy stable TBP 
thermal and radwrytic degradation in the process of uranium 
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extn^tofbnndi^pbospn^ 

fenns strong extractable complexes with ptatonium and aromhiin, whilst mono- 
butyl plic«nbate tends to to Both products lead to complications in 
the extraction process and am ren^ Therefore 
the purity of TOP must be accurately monitored to determine the point at which 



As a further example, the method of the presem inventio 

viscosity and the viscosity values may be used to determine the percentage of 
odourless kerosene (OK) in TOP. The measured viscosity is largely independent of 
the crystal used and of the droplet size. Four cxiHwnents were condiicted, each with 
a different crystal and with no particular control of droplet size. The -actuaT 
viscosity was also measured using a reverse flow viscometer in a constant 
temperature water butt 

The results are shown in Figure 7. A clear difference in measured viscosity is 
obiamedasthepenxntagetfOK Measurement of viscosity bytheinethod 

tfthepreseminvenikmcante 

whether this TBP/OK solvent system is going "c^ of specification" during its use in. 
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CLAIMS 

1. A method for determining a physical or chemical characteristic of a liquid, the 
method comprising: 

S depositing said liquid on a surface of a crystal of a quartz crystal 

microbalance to form a droplet on said surface; 
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30 



measuring a signal output of the crystal microbalance; and 

;to< 



2 a method according to Claim 1 wherein the liquid characteristic is the 



3. A method according to Claim 1 wherein the liquid characteristic is the 
IS viscosity of said liquid. 

4. A method according to any of the preceding churns, wherein said liquid is 
comprised of w single unknown chemical or a mixture of unknown chem ical s 



5., A method according to any of the preceding claims wherein the liquid 



20 



6. A method according to any of the preceding claims, wherein said droplet has 
a volume of from 0.5 to 1 ul. 

25 

7. A method according to Claim 6, wherein the volume of the droplet is about 
lul. 



8. A method according to any one of the preceding claims, wherein the signal 
output is one dependant on the evaporation rate of liquid from the crystal 
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9. A method according to any one of the preceding claims, wherein tbe signal 
output gives information on at least one physical property of the crystal 
microbalance. 

5 

10. A method according to Claim 9, wherein the at least one physical property is 
selected from the resonant frequency, the admittance, the quality factor, the 
radiofrequency voltage and the phase. 

10 11. A method according to any of the preceding claims, wherein the signal output 
is analysed by comparing values derived from the signal output with known 
values contained in a t 



15 



12. A method according to Claim 11, wherein die values derived from the 
measured values arc obtained from equation (3) as herein set out and defined. 

13. A method according to any of die preceding claims, wherein the rate of 
e vapora tion is con t rollable. 



20' 14, A method according to any of tbe preceding chums, wherein the rate of 

i is controlled by varying the te mp er a t u re at the crystal surface. 



IS. A method according to any of the preceding claims, wherein tbe rate of 
evaporation is controlled by varying tbe pressure at the crystal surface. 



16. An apparatus for determining the chemical composition of a liquid, tbe 
ap p ara tus compri sing; 
a quartz crystal microbalance; 

means for depositing said liquid on the surface of a crystal of said crystal 
30 microbalance to form a droplet on said sin lace. 

means for ca us in g the evaporation of said liquid from said surface; 
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means for measuring tbe signal output of the crystal miciobalance; and 
means for analy sin g die signal output to determine a fluid c hara ct e ristic of 
said liquid. 
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Plot of y -JniH&AtmJ fl-*xp-&) against •tapaed 
time 41 for methanol (a), ethanol (b) f propan-2-ol (c) v 
butan-4-ol (d), and peman-l-oJ (e). 

Fig. 6 
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